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Abstract. Superconductors exhibit unconventional electronic and magnetic prop- 
erties if the Cooper pair wave function breaks additional symmetries of the normal 
phase. Rotational symmetries in spin- and orbital spaces, as well as discrete sym- 
metries such as space and time inversion, may be spontaneously broken. When this 
occurs in conjunction with broken global U(l) gauge symmetry, new physical phenom- 
ena are exhibited below the superconducting transition that are characteristic of the 
broken symmetries of the pair condensate. This is particularly true of vortices and 
related defects. Superconductors with a multi-component order parameter exhibit a 
variety of different vortex structures and closely related defects that are not possible in 
condensates belonging to a one-dimensional representation. In this article we discuss 
the structure of vortices in Fermionic superfluids and superconductors which break 
chiral symmetry, i.e. combined broken time-inversion and 2D parity. In particular, we 
consider the structure of vortices and defects that might be realized in thin films of 
3 Hc-A and the layered superconductor Sr 2 Ru04, and identify some of the characteristic 
signatures of broken chiral symmetry that should be revealed by these defects. 
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1. Introduction 



The BCS theory of superconductivity pQ, combined with refinements over several 
decades [21 El SI EJ , ranks among the major achievements in theoretical physics during 
the last century. The central feature of BCS theory is Cooper pair condensation - i.e. 
the macroscopic occupation of a single quantum state of bound pairs of fermions. In 
its simplest form for s-wave, spin-singlet pairing the "order parameter", or Cooper pair 
amplitude, is given a complex scalar function of the relative coordinate, r, and center 
of mass position, R, of the pairs, 

* = (V> T (R + r/2)^(R-r/2)> = |*| e «. (1) 

The Cooper pair amplitude is both a measure of the condensate density, |\1/| 2 ~ 0(N/V), 
and a reflection of the spontaneously broken symmetry of the ordered phase, which in 
this case is the global U(l) gauge symmetry generated by Fermion particle number 
[6]. For homogeneous states the phase, is the signature of the broken U(l) 
symmetry, i.e. a macroscopic fraction of Fermions condense into a two-particle state 
with the same quantum phase. The broken global U(l) symmetry also implies a 
family of degenerate ground states related to one another by the elements of the 
symmetry group. Long-wavelength spatial variations of the phase form a branch of 
low-lying collective excitations, the Anderson-Bogoliubov phase mode [8], with an 
acoustic dispersion relation. This is a classic example of a Goldstone mode that 
accompanies a broken continuous symmetry. In superconductors, Cooper pairs are 
charged particles with charge 2e, and the requirement of local U(l) gauge symmetry 
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leads to a massive Goldstone boson via the Higgs mechanism, that corresponds to 
the longitudinal oscillation of the electromagnetic vector potential, and is a high 
energy plasmon excitation^. In neutral superfluids the phase mode is observable as 
fourth sound, a collective excitation of the condensate which has a velocity that is 
determined by the "phase stiffness", p s oc |\&oq| 2 , or superfluid density, where \1/ C q is 
the equilibrium condensate amplitude. This rigidity stabilizes the condensate against 
phase fluctuations, and is responsible for the Josephson effect j9], persistent currents 
[TO] , and flux quantization in superconductors. Indeed the standard signature of BCS 
condensation is the quantization of flux in units of $o — hc/2e, or the quantization of 
circulation in units of kq = h/2m in a neutral pair condensate of Fermions of mass m. 

Quantized flux and quantized circulation are consequences of a fundamental 
constraint on the phase of the pair amplitude, 



where N c = ± 1, ±2, ... is the winding number of the phase around a closed contour C 
within the condensate. This quantization condition reflects the physical requirement 
that the order parameter be a single-valued field, which for a complex scalar pair 
amplitude is simply that the phase return to its value modulo 2tt. For a multiply 
connected geometry this constraint and the phase stiffness implies energy barriers 
separating states of different winding number, and thus to quantized persistent currents. 
In a simply connected geometry states with iV c ^ force the condensate to effectively 
become multiply connected as there is necessarily one or more phase singularities 
interior to the contour C. These singularities, points in two dimensions or lines in three 
dimensions, form a spectrum of topologically stable "defects" of the order parameter, 
i.e. phase vortices, labelled by their winding number. These topologically stable defects 
are often energetically stable, or metastable because of energy barriers separating states 
with different winding numbers. The local structure of these defects, the vortex cores, 
are determined in part by their topology - their winding number in this simple case, the 
local spectrum of fermionic excitations and their interactions, as well as scattering of 
the Fermions by impurities or other sources of disorder [TTj . 

BCS superconductors exhibit unconventional properties if the order parameter 
breaks additional symmetries of the normal phase. Spin- and orbital rotational 
invariance, as well as space and time inversion, may be spontaneously broken. When this 
occurs in conjunction with broken U(l) gauge symmetry, new collective excitations of 
the pair condensate [12], [13] , novel heat and transport properties (THUS], unconventional 
vortices [HI |T7] as well as point defects [TH [19] characteristic of the complex symmetry 
breaking are possible [20] . In this article we discuss the structure of vortices in superfluid 
3 He and superconductors which break spin and orbital rotation symmetries, parity 
and/or time- inversion symmetry. This class of pairing states is believed to describe the 

| A direct consequence of the Higgs mechanism is the Meissner effect in superconductors. 




(2) 
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ground states of thin films of 3 He-A[2T], the low-temperature superconducting phase of 
UPt 3 [22J, and the layered superconductor Sr 2 Ru0 4 [23]. 

2. Theoretical Formalism 

The superconducting order parameter can be defined in terms of the condensate 
amplitude for Cooper pairs, 

V a p= (Mti)Mt2)), (3) 

where ip a ( r ) is the fermion field operator for spin projection a. This amplitude is directly 
related to the anomalous Matsubara propagator [24J, 

&*/3(ri,Ti,r2,T 2 ) = - (T T V>a(ri,Ti)V>0(r2,T 2 )) 

= - (T T ^(R + r/2, r)^(R - r/2, 0)) = £ Q/3 (R, r; r) , (4) 

expressed here in terms of the center-of-mass, R, and relative (orbital) coordinate, r, 
of the pair, and the (imaginary) time difference, r = n — r 2 . T T is the (imaginary) 
time-ordering operation for Fermions. Fourier transforming with respect to the orbital 
coordinate and time gives the anomalous propagator for pairs with orbital momentum 
p, center of mass position R, 

^(p,R;e n ) = J dre^ h j dr e-^&^R, r; r) , (5) 

and Matsubara energy, e n = {2n + l)vr/ /5 for n = 0, ±1, . . ., where j3 = l/ksT. The 
local spectral function for the pairs is obtained by analytic continuation to the real axis, 
ie n -^e + i0 + [21]. 

The orbital radius of a Cooper pair is of order £ = hvf /27t/cbT c , which is typically 
large compared with atomic scales, e.g. the Fermi wavelength, £o 3> fi/pf- In this limit 
we can factor out the short wavelength variations of the pair amplitude and compute 
directly the spatial variations of the pair amplitude and order parameter on mesoscopic 
length scales defined by the coherence length, £ 0; the mean-free path, I and/or the 
London penetration depth, A. This separation occurs because the high-energy, e ~ Ef, 
short-wavelength, |p| ~ pf, properties of the Fermi liquid are unaffected by the long- 
wavelength, low-energy pairing correlations to leading order in h/pf^ , h/pji, ksTc/ Ef 
etc. Thus, we define the amplitude for Cooper pairs near the Fermi surface by integrating 
over the low-energy band of states near the Fermi level [26] , 

d£p&*/9(p,R;en), (6) 

-n c 

where pj is the Fermi momentum, £ p = f/(|p| — Pf) and Q c Ef is the quasiparticle 
bandwidth. A detailed discussion of this procedure is given in Refs. [25], [261 E3 EHJ [29] . A 
synopsis of the quasiclassical transport theory for equilibrium states of inhomogeneous 
superconductors, as well as the methods used to calculate the order parameter and 
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electronic structure of superconducting states associated with impurities, interfaces and 
vortices is described below. 

The propagator for Cooper pairs with relative momenta on the Fermi surface 
given in Eq. [6] is coupled to the low-energy propagator for fermionic quasiparticles 
of the superconductor. The coupled equations for quasiparticles and Cooper pairs are 
formulated in terms of a 4 x 4 matrix propagator in the combined particle-hole (Nambu) 
and spin space, which incorporates the the spin correlations that are required by the 
pairing correlations, 

S(p, R; e n ) =-J Q P dT e fc » T J d 3 r e^ (T^n, r)^(r 2 , 0)) , (7) 

where $ = (ip^ , ip^ , ipt , ip}) is the four-component (Nambu spinor) Fermion field 
operator, i"i = R + r/2, r 2 = R — r/2. Note also that for imaginary time evolution, 
-0(r, r) = ip T(r, — r). The quasiclassical matrix propagator is then defined by integrating 
the over the low-energy band of states near the Fermi surface, |£ p | < Q c , with 
k B T c < tt c < Ef, 

$(p/,R;e n ) = - ! + C d£ p ? 3 %(p,R;e n ) = (f [ ) , (8) 

where t% is the third Pauli matrix in Nambu space, and we have renormalized by dividing 
by the spectral weight, a, of the normal-state quasiparticle resonance, < a < 1. Its 
structure in Nambu space is represented by the quasiparticle propagators, q and g, and 
the pair propagators, f and f. Each of these propagators is a 2 x 2 matrix in spin 
space denoted by small hats. The pair propagator naturally separates into spin-singlet 
(anti-symmetric) and spin-triplet (symmetric) amplitudes, 

fap(Pf, R; ^) = fo(p/, R; e«) (^) a/3 + f(p/, R; <*) ■ (ivo- y ) af3 (9) 

while the natural description for the quasiparticle propagator is in terms of spin-scalar 
and vector components, 

a/3 (p/, R; e n ) = g(p/, R; e n ) 5 a/3 + g(p /; R; e n ) ■ a aP , (10) 

where a = {o~ X) a y) a z ) are the Pauli spin matrices. The description in terms for four 
matrix propagators is convenient, but redundant. Fundamental symmetry relations 
associated with Fermion anti-symmetry and conjugation (i.e particle <-> hole) connect 
the two diagonal and off-diagonal propagators, 

s(p/, R; e n ) = s(-p/, R; e„)* = g(-p /; R; -e n ) tr (11) 

?( P/ , R; e n ) = f(- P/ , R; e n f = -f(p/, R; -e„) f . (12) 

where b tr (b^) is the transpose (adjoint) of b. 
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2.1. Transport Equation 

Extensions of Landau's theory of normal Fermi liquids to include BCS pairing 
correlations [30], electron-phonon interactions [31], and disorder [32] culminated in the 
equations of Eilenberger [5] and Larkin and Ovchinnikov [33J , formulated in terms of the 
Nambu matrix propagator, g(p/,R;e„), obeying transport-type differential equations 
along classical trajectories defined by the Fermi momentum J§| 



where [. . .] is the commutator in Nambu space between the matrix propagator, g, 
the (imaginary) time-development operator, %d T <-> ie n %, and various internal and 
external fields that couple to quasiparticles and pairs. Eilenberger's transport equation 
is obtained by an energy- and momentum-space renormalization of Dyson's equation 
and the full many-body equations for the self-energy functional. The normalization 
information contained in Dyson's equation is absent in Eilenberger's equation, but is 
recovered in the form of a contraint [5], [33] that must be satisfied by the physically allowed 
solutions of the transport equation. The normalization constraint on the Matsubara 
propagator is 



This constraint can be obtained as a boundary condition which supplements the 
transport equation, and reflects the fact that in the absence of any disturbance the 
propagator reduces to that for local homogeneous equilbrium. For more detailed 
discussions of the normalization condition see Refs. [5l |33| [36]. 

The renormalization procedure leads to an expansion of the full many-body Green's 
functions and self-energies in terms of ratios of low to high energy scales, e.g. ksTc/Ef, 
u cxt /Ef, or short to long wavelength scales, e.g. VP/6o> ^l/Pf- We use a single 
parameter, small, to classify the order of magnitude of various terms in the transport 
equation, the magnitude of the self-energy and external fields, etc. The results reported 
here are based on the leading order expansion of the self energy in small, particularly the 
mean-field pairing self-energy, A, that describes the inhomogeneous equilibrium state of 
the superconductor, and the effects disorder in terms of scattering by a dilute random 
distribution of impurities represented by the impurity self energy, S. However, we omit 
the Landau molecular field self-energy, which is also of order small. We do not expect 
the Landau molecular field to lead us to qualitatively different conclusions regarding 
the basic structure of vortices in chiral p-wave states, however these terms may play a 
role in the relative stability of various multi-vortex configurations. We also treat the 
mean-field pairing self energy in the weak-coupling limit. Indeed a useful operational 

§ Extensions of the transport theory of Eilenberger to non-equilibrium superconductivity were 
formulated by Eliashberg |34j and by Larkin and Ovchinnikov [35] - The nonequilibrium transport 
theory is discussed our companion paper on vortex dynamics in this volume [40j . 



e n f 3 — A — E , g = , 



(13) 




(14) 
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definition for the order parameter is the mean field pairing self energy, 



where A is the 2x2 spin matrix order parameter with components A a/ j(p/, R; e n ), 
loosely referred to as the "gap function", is is both a measure of the pairing correlation 
energy and the broken symmetry of the superconducting state. The gap function is 
determined by the BCS self-consistency equation relating A(pj,R;e n ) to the Cooper 
pair amplitude, f(p/,R;e n ), and the interaction, A(p/, e n , p'f, e' n ), that provides the 
"pairing glue". In the weak-coupling limit the frequency dependence of the pairing 
interaction, and the pairing self-energy, is constant within the bandwidth |e n | < Q c . 
The order parameter then satisfies the weak-coupling gap equation, 

A Q/3 ( P/ R)= P f ,a Pf ,(3 

, Kl<n c 

= N f / d 2 p' f x a ^ p (p f , P ' f )T f7p(p/;0 ( 16 ) 

J e' 

where f = — «-» — is the Cooper pair propagator and is the bosonic propagator 

responsible for the pairing glue. Nf is the single-spin normal-state density of states at 
the Fermi level and the integration is an average over the Fermi surface normalized to 
J d 2 pj = f dS Pf n(pf) = 1, where n(pf) is the normalized angle-resolved density of 
states on the Fermi surface. 

The impurity self-energy is a functional of the quasiclassical propagator and defined 
in terms of the self-consistent impurity scattering t-matrix, 

S(p/,R;e n ) =n s t(p / ,p / ,e„;R) , (17) 

where n s is the density of impurity density and t(p/, P/, e„; R) is the forward scatting 
limit of the self-consistent t-matrix describing multiple-scattering of quasiparticles and 
pairs by an impurity in the superconductor, 



t(p/,P/;e n ) 



+ 



= u( P/ , p' f ; e n ) +N f J d 2 p" f u( P/ , pj; e n ) fl(pj; e n )t(pj, p' f ; e n ) (18) 

where fl(p/,e n ) = = is the full matrix propagator and u(pj,p'j), represented by 
the cross and dotted line, is the matrix element of the impurity potential between 
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normal state quasiparticles with momenta p/ and p'^ on the Fermi surface. We 
model the disorder in terms of non-magnetic impurities that scatter quasiparticles 
isotropically, i.e. we retain only the s-wave contribution to the impurity potential, thus, 
u = uq%. The corresponding scattering phase shift in the s-wave channel is defined by 
So = tan^ 1 (77 NfU ). The t-matrix equation, and thus the impurity self-energy, can be 
expressed as the solution of the matrix equation, 

t(e n ) =u + N / u(5(pJ;e n )>t(e n ), (19) 

where (jj(p^;e n )) is the Fermi-surface average of the propagator. In the normal state, 
the low-energy quasiclassical propagator reduces to q n = —iir sgn(e n )%, and thus the 
t-matrix descbribing the scattering of normal state quasiparticles and quasiholes near 
the Fermi surface becomes, 

t N (e n ) = 4r sin 5 ° e- iS S n ^ Sof * . (20) 

In this model the disorder is characterized by the strength of impurity potential, uq, the 
mean density of scatterers, n s and the density of states for quasiparticles at the Fermi 
surface. The mean lifetime of a quasiparticle on the Fermi surface in the normal state, 
r N , is obtained by analytic continuation of the self-energy to the real axis to obtain the 
retarded self-energy, T^(e) = S N (e n — > e + i0 + ), 

-**™ = wr (21) 

The mean free path, £ N = VfT N , reduces to the classic result, £ N = \jn s a, where a is the 
total cross section for scattering of a quasiparticle off an impurity, a = (4irh 2 /p'j) sin5 2 . 
In terms of the dimensionless cross section a = sin 2 Sq, the weak scattering limit ("Born 
limit") corresponds to a — > 0, while the strong scattering ("unitarity limit") corresponds 
to a — 1. 

For any inhomogeneous superconducting state, the relative importance of impurity 
scattering is determined by two parameters, (i) a pair-breaking parameter defined by 
ratio of the scale of the pair correlation length and the mean free path, x = (,q/£ n and (ii) 
the strength of the impurity scattering potential as measured by the dimensionless cross 
section, a. For unconventional pairing states, such as chiral p-wave superconductors, 
disorder destroys the superconducting state at a critical level of disorder given by 
x c ~ 0.28. The results reported below are based on self-consistent calculations of the 
order parameter, current distribution and local density of states for the clean limit and 
for modest levels of disorder with £ N = 10£o in the Born limit. 

2.2. Observables 

To complete the theoretical framework for our calculations of the vortex structure for 
Fermi-liquid superconductors and superfluids we include expressions for the equilibrium 
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current, and local density of states in terms of the propagator. In particular, we define 
the normalized angle-resolved local density of states (LDOS) by dividing out the the 
angle-resolved normal-state DOS at the Fermi level, iV/n(p/). Thus, the normalized 
local spectral density due to particle-hole coherence of the pair condensate and the 
formation of sub-gap excitations associated with scattering and spatial inhomogeneities 
of the order parameter is obtained from the retarded propagators for particle and hole 
excitations by analytic continuation to the real axis, 

JV( P/ , R; e) = -i-JmTr{? 3 R (P/, R; e)} , (22) 

where the trace is over the 4x4 Nambu space. All equilibrium properties are determined 
by this spectrum. In particular, the equilibrium charge current associated with a 
topological defect in a superconductor is given by, 

/r+oo 
d 2 p f / de [ev f (pf)] Af(p,, R; e) (2/(e) - 1) (23) 
J — oo 

= N f f d 2 p f [ev/fcv)] ^Tr{f 3 g( P/ ,R;e n } , (24) 

where /(e) = l/(e /3e + 1) is the Fermi distribution. The orbital magnetization density 
is then given by V x m = — - j(R). For equilibrium states it is convenient to calculate 
the total current by transforming to the Matsubara representation, as shown in Eq. [241 
However, if we are interested in the contributions to the current from specific regions 
of the excitation spectrum then the formulation in terms of the LDOS is more useful. 
In particular, the angle-resolved spectral current density is defined as the net current 
density carried by states at ±p/ and energy e [TTj . 

j( P/ , R; e) = N f [ev f ( P/ )] [jV( P/ , R; e) - JV(- P/ , R; e)] ■ (25) 



Ricatti Equations 

The order parameter and self-energy must be determined self-consistently with the 
solution of the transport equation for the propagator. An efficient method for solving the 
transport equation is based on a parameterization for the propagator that automatically 
satisfies the normalization constraint, and thus eliminates the possibility of spurious 
solutions. The parameterization is defined by [37J EH USE IZE] 



-27 




g = -ivr^rt I 1 + T T (26) 



where the prefactor is given by 



(1-77)- 1 
(1-77)" 
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The coherence amplitudes 7 and 7 are 2x2 matrices in spin space which obey Ricatti- 
type equations, 

iv f ■ V7 + 2ie n <7 - 7A7 - 2X7 + A = , (28) 
iv f ■ V7 - 2ie„7 - 7A7 + 2S7 + A = , (29) 

and are simply related to the particle- and hole-like projections of the off-diagonal 
propagator, 

7 = -(i 7 r-fl)- 1 f, 7 = (m + i)- 1 ?, (30) 

where the projection operators in Nambu space for particle-like (+) and hole-like (-) 
excitations are given by [36] . 




In particular, = P± and P+P- = P_P + = follow directly from the normalization 
condition. The physical interpretation of these projectors, and the Ricatti amplitudes, 
is easily established by comparing the action of the projection operators on a general 
Nambu spinor with the particle-like and hole-like solutions of the Bogoliubov or Andreev 
equations. The fundamental symmetry relation in Eqs. [1X11121 for the particle and hole 
components of the quasiclassical propagator also imply symmetry relations relating the 
two Ricatti amplitudes, 

7(p/,R;e n ) = 7 (- P/ ,R;e n )^. (32) 

The utility of the Ricatti equations is that they provide an efficient approach to solving 
the transport equation. The Ricatti equations are easily integrated numerically because 
they have numerically stable solutions. The extension of this formalism applicable to 
vortex dynamics is described in a companion article for this special collection, Ref. |40j, 
as well as Refs. [381 [391 E] . 



2.3. Calculational Methods 

The calculations for topological defects reported here are obtained from self-consistent 
solutions for the propagator, order parameter and self energy. For two-dimensional 
layered superconductors or superfluids we define all relevant functions, e.g. g(p/, R, e n ), 
on a discrete set of points in both position and momentum space. The center of mass 
coordiate is discretized as R = h(ma + mb), where h is the lattice spacing, a and 
b are lattice unit vectors and (m, n) is a pair of integers defining a discrete lattice 
point. A typical choice for the lattice spacing is h = 0.25 £0, and the dimensions of 
the largest grids used to carry out these calculations was 120 x 120 lattice sites. The 
larger grids were used to investigate the stability of vortices multiple units of phase 
winding. We also used both rectangular and hexagonal lattices. For calculations of 
isolated defects we place the defect initially at the origin and enforce the constraint 
on the global phase winding by fixing the order parameter by its asymptotic form for 



Vortices in chiral, spin-triplet superconductors and superfluids 



11 



|R| — > oo along an exterior contour as shown in Fig. [TJ At each grid point we use a 
fourth order Runge-Kutta algorithm to integrate the transport equation along a discrete 
set of trajectories representative of the Fermi surface and the basis functions defining 
the pairing symmetry. For the calculations reported here we used sixteen trajectories 
equally spaced in momentum space around the Fermi surface, which for simplicity we 
assumed to be cylindrical. However it is straight forward to implement more detailed 
Fermi surface geometries if needed (see Ref. [12] )• We also calculate the propagator 
for a finite set of Matsubara frequencies. The maximum Matsubara frequency required 
to achieve high precision increases with decreasing temperature. For T/T c = 0.4 eight 
Matsubara frequencies is typically sufficient. 



j> V-d ■ d£= m x 2tt 
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Figure 1. Computational grid with spacing h. At each grid point, labelled by 
R = h(m, n), we integrate the transport equation along a classical trajectory (magenta) 
defined by the Fermi velocity, V/. The global phase winding is enforced on the outer 
contour. The defect (red dot) is initially positioned at the center of the grid. 

A calculation is initialized by a "seed" for the order parameter field, with the 
topological constraint implemented on the boundary as described above. The propagator 
is then calculated for all the grid points, and for the set of Fermi surface trajectories 
and Matusbara frequencies at each grid point. The order parameter and impurity self 
energy are then calculated from discretized forms of Eqs. [JEl [UJ and [191 The updated 
values of the order parameter and impurity self-energy are then used as new inputs to 
the transport equation, which is solved again for each grid point to obtain an improved 
solution for the propagator. This iterative procedure is continued until it converges to 
a specified precision for the deviations of the order parameter and self energy between 
iterations. The final results for the equilibrium order parameter and excitation spectrum 



Vortices in chiral, spin-triplet superconductors and superfluids 



12 



are found to be insensitive to the initial order parameter field, except of course for the 
constraint on the phase winding. 

2.4- Pairing Symmetry 

Fermi statistics requires that pair amplitude, f ai g, the pairing interaction, A Q/ g ;7P , and 
thus the gap function, A a/ g, obey the anti-symmetry condition jjj] 

Aa / j(p / ,R) = -A /9a (-p / ,R). (33) 

For materials in which the normal metallic state possesses inversion symmetry 
the pairing interaction separates into even- (g) and odd-parity (u) channels, which are 
respectively anti-symmetric and symmetric under spin exchange and correspond to the 
spin-singlet (S = 0) and spin-triplet (S = 1) pairing channels, 

<-*■(") 

Aa/3 ;7P (p/, p' f ) = {i^ y )a(3 A (9) (P/> p'f) {iv v ) 1P + fay^a/i • A (p/, p' f ) ■ {iaa y ) ip , (34) 

Furthermore, the pairing interaction separates into a sum over invariant bilinear 
products of basis functions for the irreducible representations, T, of the point group, 
for both even- (f/rv(p/)) and odd-parity (?yrv(p/)) sectors, 

A (S) (P/>P/) = E A ^,(P/)C(P/), (35) 

Vg,V 

«(«) 

A (P/, p' f ) = Ar^r,(p/) ® C(P/) • ( 36 ) 

These basis functions are the eigenfunctions of the linearized gap equation, while A r is 
the eigenvalue that determines the instability temperature for Cooper pairs belonging 
to the the irreducible representation, T, i.e. 1/Ar = In(1.13f2 c /T c r ). Except in rare cases 
of accidental degeneracy, or a weakly broken symmetry [16], the pairing interactions 
are well separated and the most attractive channel determines both the transition 
temperature, T c and the basis functions, {r] rv }, for the irreducible representation, T, 
that determines the pairing state(s) of the superconductor. The index v labels the 
basis functions for the representation T with dimension d r . Although the exact basis 
functions depend on details of the anisotropic Fermi surface and pairing interaction, 
representative basis functions which exhibit all the broken symmetry properties of the 
irreducible representation are easily constructed [17]. Below we list the irreducible 
representations and representative basis functions for the group D^, appropriate for 
Sr 2 Ru0 4 in Table El 

Unconventional pairing states which have no equal-time average, i.e. 'odd- frequency' pairing states, 
obey a more general anti-symmetry condition which includes the imaginary time coordinate [43] . Such 
states have the parity assignments for total spin quantum numbers interchanged. See Berezinskii 44 , 
and for a more recent discussion of these states, Balatsky and Abrahams [45] , 
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Table 1. Representative basis functions for D^h appropriate to strong spin-orbit 
coupling with d||c, or negligible spin-orbit coupling and arbitrary direction of d. 



Thus, barring accidental near degeneracy [48], the order parameter, A a/3 (p), is 
defined by a single representation, and is either even- or odd-parity, and therefore spin- 
singlet or spin-triplet, respectively{p| 

A Q/3 (p/) = A(p/) (io- y ) a /3 , singlet (5 = 0), (37) 

A a/3(P/) = A(p/) • (iao- y ) a p , triplet (5=1), (38) 

with A(pj) = A(— pf) and A(p^) = — A(— pj). The expansion in the eigenf unctions 
belonging to the dominant representation T (either even or odd-parity) of the linearized 
gap equation gives, 

A( P/ )= ^A^(P/), (39) 

V 

d T 

A( P/ )= ^A^( P/ ). (40) 

v 

2.5. Odd-Parity, Spin-Triplet Pairing 

The prototype for odd-parity, spin-triplet pairing is realized in the neutral Fermi 
superfluid 3 He. In particular the A-phase of superfluid 3 He is identified as an equal-spin 
pairing state with a p-wave order parameter of the form 



A( P/ ) = ^rf(m + m)-p, (41) 

where i = m x n is the quantization axis for the orbital angular momentum of the 
pairs, i.e. £■ L° rb (m + in) ■ p = +h (m + in) ■ p, p = p//|pj|and d is the axis along 
which the pairs have zero spin projection, i.e. d ■ S pair A = 0. This corresponds 

P For materials in which spin-orbit effects are strong, i.e. actinide and rare earth heavy fermion metals, 
the labels " spin-singlet" and " spin-triplet" do not refer simply to the eigenvalues of the spin operator 
for electrons, but rather a "pseudo-spin" . The Kramers' degeneracy in zero-field guarantees that each 
p state is two-fold degenerate. These two states may be labeled by a pseudo-spin quantum number a. 
Thus, for many of our considerations the distinction between spin and pseudo-spin is not important 

\MMM- 
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to a spin-state that is an equal amplitude superposition of | ft) an d | ||) spin 
states in the plane perpendicular to d. In bulk 3 He-A, which is stable over a narrow 
temperature range at high pressures, both the orbital and spin quantization axes are 
broken symmetry directions, aligned only by external walls or weak symmetry breaking 
fields. In particular, the very weak nuclear dipolar energy is minimized by orienting 
the nuclear spins in the orbital plane, i.e. d\ \ ± £ [53]. In thin films of superfluid 3 He 
scattering by the surface and substrate leads to substantial pair-breaking unless the 
pairs are confined to the x-y plane of the film [51]. This condition favors the A-phase 
with <i||£||z over the B-phase for the entire pressure range (c.f. Fig. 3 in Ref. |21j). 
However, in this reduced geometry the planar phase, which is also an equal-spin pairing 
state, but with d z = is degenerate with the A-phase in weak-coupling theory. It is 
unknown whether or not the A-phase or the planar phase is the stable phase in thin 
films, particularly at low pressures where strong-coupling effects are extremely small 
[55] . For the planar state the two spin states are aligned out of the plane of film 
and correlated with the orbital pairing state to form a state with zero total angular 
momentum projection along z, 

A 

A(P/) = — {(x + iy) (p a - ip y ) + (x- iy) (p x + ip y )} . (42) 
2.6. Broken Time-Inversion and Chirality 

In contrast to the planar phase the A-phase spontaneously breaks both time-reversal 
symmetry and reflection symmetry in a plane normal to the film. This combination 
of broken 2D reflection symmetry and time-reversal symmetry is referred to as broken 
chiral symmetry. The structure of vortices in chiral superfluid 3 He-A are fundamentally 
different than vortices in non-chiral condensates such as the planar phase of superfluid 
3 He. 

There are also candidates for triplet pairing in the class of strongly correlated 
electronic superconductors. For example, UPt3 is believed to have an odd-parity order 
parameter belonging to the E 2u representation of D 6h with a low-temperature phase that 
is an equal-spin pairing state with d\\z, and breaks parity and time-reversal symmetry, 
i.e. A = dp z (p x + ip^) 2 The evidence and analysis favoring this identification comes 
from observations of a multiple superconducting phases, low-temperature transport 
properties and observation of anisotropic Pauli limiting, and is summarized in Refs. 
[56| |22~1 ITS] . More recent small-angle neutron scattering studies of the flux lattice 
transitions also support this identification [57| 158]. 

The layered superconductor, Sr 2 Ru0 4 is also candidate for chiral spin-triplet pairing 
with a proposed order parameter that is essentially identical to that for a thin film of 
3 He-A [59]. Strong evidence for unconventional pairing in Sr2Ru04 is provided by 

+ The higher symmetry of the hexagonal rotation group allows for a second two-dimensional E 
representation formed from the basis functions of the two B representations of D^. The odd-parity 
version, E2„, with d\\z has basis functions dp z (p x ± ipy) 2 ■ 
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the suppression of T c below 1.5 K by non-magnetic impurities [60J, while the absence 
of a Hebel-Slichter peak in \jT{T probed by NQR and Knight shift measurements 
[61] are cited in support a spin-triplet order parameter with the spin-quantization axis 
d\\z. Phase-sensitive measurements based on the current-phase relation for Josephson 
junctions coupling a conventional s-wave superconductor, Au-In alloy, and Sr 2 Ru04 
support the interpretation that Sr 2 Ru0 4 is an odd-parity superconductor [62]. Studies 
of the excess current obtained from point contact spectroscopy agree quantitatively 
with a p-wave triplet pairing state [63]. Evidence for broken time- reversal symmetry 
of the superconducting order parameter is provided by yuSR [M] and magneto-optical 
Kerr effect [65J experiments indicating that spontaneous currents and magnetic fields 
develop below T c in the Meissner state. However, the identification of broken chirality 
is controversial because local magnetic probes have not observed the fields associated 
with the domain- wall currents that are expected from broken chiral symmetry [66J. 

3. Vortex structure in chiral p-wave, spin-triplet condensates 

Superconductors with a multi-component order parameter exhibit a variety of 
different vortex structures not possible in condensates belonging to a one-dimensional 
representation. These structures can lead to multiple superconducting phases as well 
as novel electronic and magnetic properties [17J. Theoretical studies of these vortex 
structures can provide important information for future experimental studies of vortex 
phases in triplet superconductors. 

Theoretical analysis of the stability of the topologically stable vortex structures 
requires calculations of the free energy of these structures, typically taking into account 
the interactions with other vortices in the vortex lattice structure. Here we restrict our 
analysis and discussion to the internal structure of individual vortices appropriate to 
the low- field regime H ~ H C1 . 

For any inhomogeneous superconducting state, and specifically for a vortex in a 
p-wave, spin-triplet superconductor or superfluid, the order parameter can be expanded 
in basis functions of the relevant irreducible representation, T, of the symmetry group, 

A(p /l R) = 53A v (R)^ i „(p / ). (43) 

For chiral p-wave states appropriate for thin films of 3 He-A and Sr2Ru04 with d\\z, 
the representation is two-dimensional and the basis functions, ffr, u {Pf) — dr}„(pf) with 
T)±{pf) = p x ±ip y , correspond to orbital angular momentum eigenstates with L°* h = ±h. 
For 3 He-A these basis functions are exact to the extent that higher-order representations, 
e.g. f-waves with L z = ±1, can be negelected [67J. In the case of S^RuC^ higher-order 
harmonics belonging to the same irreducible representation, E u , are possible. Typically, 
we can neglect these higher-order terms as the basic features of the vortex structure are 
obtained with the lowest order odd-parity harmonics. 

Because time-reversal symmetry is broken for chiral, p-wave superconductors, there 
are two possible degenerate ground states defined by their spontaneous orbital motion. 
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The homogeneous ground state is spontaneously chosen between the two degenerate 
time-reversed states with L° rb = ±h [16]. In the following sections we consider vortices 
in a single domain with internal orbital angular momentum L° rb = +h. Vortices in 
the time-reversed ground state can be obtained simply by applying the time-inversion 
operation: A — * T * A. 



3.1. Global and local winding numbers for vortices 

Consider an isolated vortex with circulation m x 2ir far from the vortex core, with m 
an integer. In this limit the amplitude approaches its bulk value but carries the global 
phase winding, 

lim A( P/ , R) = A^pf, R) = d ^ e im4 > (p x + ip y ) . (44) 

|R|->oo \/2 

Asymptotically the vortex is cylindrically symmetric, and is thus an eigenfunction of 
the generator, L z , for rotations about the z axisQ 

L z * Aoo(p/, R) = / ft Aoo(p/, R) , Z = 0,±l,±2,... (45) 

where L z = L™ + L° rb is the sum of the z-component of the angular momentum operator 
for the center of mass (CM) of the pairs, = |<9/<90 with R = |R|(cos0, sin0), and 
the orbital angular momentum of the pairs relative to the CM, L c z m = ^d/dip p , where (p p 
is the azimuthal angle of the relative momentum of the pairs, p/ = pj(cosLp p , cos(p p ). 
Thus, we have I — m + 1. At finite, but large distances from the vortex core, i.e. 
| R.| ^> £o, gradients of the asymptotic vortex order parameter induce the degenerate, 
time-reversed pairing state leading to an order parameter of the more general form, 



A( P/ ,R)=d 



|A i (R)l e im4> l* y > + |A_(R)| e ip * 



V2 V2 



(46) 



It is clear that the time-reversed phase with winding number p is nucleated in the 
vicinity of the vortex core. At large distances the vortex retains axial symmetry. This 
enforces a condition on the local winding number, p, of the time-reversed phase in the 
core, i.e. p — 1 —m + l [22]. This condition is exact for vortices in a pure p-wave chiral 
ground state such as 3 He-A. For chiral ground states in superconductors with discrete 
symmetry the condition on the local phase winding can be generalized. For example, for 
the discrete point group D 4/l appropriate to Sr 2 Ru0 4 the condition on the local phase 
winding becomes, 

p = m + 2 + An , (47) 

where n = 0, ±1, ±2, . . . includes the effect of higher order harmonics comprising the 
ground state that are invariant under 4-fold rotations. If deviations from cylindrical 
symmetry are small we can neglect corrections coming from n ^ 0. For simplicity we 
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Table 2. Table of vortex states for the chiral ground state with L° lh = +H. Shown are 
the dominant order parameter, A+, its global phase winding to, and the corresponding 
vortex core order parameter, A_, and its local phase winding p. The color code 
describes the phase winding - a change in color corresponds to a change of the sign in 
the real part of the corresponding order parameter amplitude. The order parameter 
with to = and p = +2 corresponds to a defect with no global phase winding, but 
with local phase winding near the defect. 

restrict our quantitative analysis to vortices with n — 0, but comment on the qualitative 
effects of the higher order harmonics where appropriate. 

Table 13.11 summarizes the lowest energy vortex states and their corresponding 
winding numbers. The figures show calculations of the dominant order parameter, 
A + , as well as the sub-dominant order parameter, A_, that develops in the core of 
the vortex. The color coding indicates the corresponding phase winding - a change in 
color corresponds to a change of the sign in the real part of the corresponding order 
parameter amplitude. The size of the circular symbols indicates the relative magnitude 
of the order parameter. We discuss each of the vortex states in more detail below. 

Figures El [3J and [6] provide additional quantitative results for the order parameter 
structure, the current distribution and local density of states for the vortices listed in 
Table [37T1 The magnitudes of the order parameter components, A + and A_, are shown 
in column 3 as a function of distance along the x axis through the vortex center. The 
local density of states, A/"(p/,R;e) is shown for a family of trajectories parallel to the 
x axis with impact parameters y separated by Ay = 0.78£o- Measurements of the local 
density of states using STM spectroscopy, which has been successfully applied to vortex 
studies in layered superconductors [68J including the high-T c cuprates [69], could provide 
valuable information on the vortex structure and pairing symmetry of S^RuO^ 

Our calculations were performed with impurity scattering included in Born limit 
for a mean free path of i = 10£o, where £ = hvf/27ikBT c is the coherence length. The 
Fermi surface parameters are assumed to be isotropic, and the temperature was chosen 
to be T = 0.2 T c . For simplicity we also assumed the high-K limit, where the penetration 
depth is large compared to the coherence length. Our calculations of the vortex structure 

* Spin rotations play no role for 3 He-A films with fixed d\\z . 
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are appropriate to the low-field limit near if cl where vortices are well separated from 
each other, or for H < if cl in the case of the defect with m = 0. The order parameter, 
impurity self energy, and the equilibrium spectra and current densities were obtained 
self consistently using the Riccati formulation of the quasiclassical transport theory with 
impurity and pairing self-energies. 

3.2. Two types of singly quantized vortices 

Broken time-reversal symmetry of the ground state leads to an interesting effect: vortices 
with opposite global phase winding numbers are inequivalent jTTj . The ground state 
pairs have internal orbital angular momentum, L° rb , which can be either parallel or 
antiparallel to the external field H = Hz that nucleates a vortex and fixes the sign of 
the global phase winding. For H||L 2 = +hz the m — l,p — 3 vortex is realized, whereas 
for the antiparallel case, — H||L 2 = +hz or H||L 2 = —hz, the m = — l,p = +1 vortex 
is nucleated. The structures of both vortices are shown in Fig. [2j Both vortices exhibit 
the large zero-energy Andreev bound state at the vortex center that is characteristic 
of odd winding number [UJ. Note however the difference in magnitude of the induced 
components for m = ±1. The relative orientation of the field, H, and the spontaneously 
chosen direction for the internal orbital angular momentum of the pairs in the ground 
state, L° rb = ±hz, determines which type of singly quantized vortex is nucleated from 
the Meissner state at the lower critical field. The m = +1 vortex has a local phase 
winding of p = +3, and thus different core energy than the simpler m = — 1 vortex with 
p — +1. As originally noted by Tokuyasu et al. [17], the difference in free energy for 
these two vortices leads to a splitting of the lower critical field, H Cl+ ^ H Cl _, for fields 
parallel (+) or antiparallel (-) to the internal orbital angular momentum of the pairs, 
L° rb [17] • Observation of this asymmetry would provide a direct signature of broken 
chirality in the ground state. 

Ginzburg-Landau (GL) calculations by Tokoyasu et al. [17] also predicted the 
possibility of spontaneously broken axial symmetry by the cores of singly quantized 
vortices with m = ±1. In the GL theory whether or not the axial symmetry is 
spontaneously broken is determined by competition between the suppression of the 
condensation energy by the global phase winding and the internal Josephson phase- 
locking energy between the two time-reversed order parameter amplitudes. If the GL 
coefficient, fa, that determines the coupling energy, F c = 4/3 2 |A + | 2 |A_| 2 , is sufficiently 
weak compared with the condensation terms, (3\ [|A + | 4 + |A_| 4 ], i.e. < /3 2 < 
then axial symmetry is spontaneously broken by the vortex core. Our results for singly 
quantized vortices, which do not break axial symmetry at any temperature, also agree 
with those of Ref. [TTJ since our theory reduces to the GL equations with the weak- 
coupling value of P2/P1 — \ i n the limit T w T c . 

Other authors have also considered differences between singly quantized vortices 
with opposite phase windings relative to the chirality of the ground state. The authors 
of Ref. [70J report results for the field dependence of the vortex structure in chiral p- 
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wave superconductors. They argue that H C2 differs for states with opposite chirality, i.e. 
L° z Ib = ±h. However, this is incorrect. The upper critical field represents a second-order 
transition between the normal state and the superconducting vortex state. Since chiral 
symmetry is unbroken in the normal state, the second-order instability field, H C2 , is 
necessarily the same for H| | + z or H| | — z. Furthermore, at H C2 both chiral components 
of the 2D representation are nucleated. As the field is reduced, the vortex density 
decreases and eventually a single chiral domain may remain, at least in the idealized 
limit of a perfect crystal with no vortex pinning. The chirality of the resulting ground 
state in the limit H — > will be determined by the direction of H, but H C2 is the same 
for either field orientation. The situation is completely different for H C1 . In this case a 
vortex enters the superconductor from the Meissner state which spontaneously breaks 
chiral symmetry. The relative sign of the phase winding of the vortex and the chirality of 
the pre-existing Meissner phase lead to different vortex core states nucleating depending 
on whether H||+zorH|| — z, and a corresponding asymmetry in H C1 . 

Recently, Yokoyama et al. [71] proposed that broken chiral symmetry should be 
observable as an asymmetry in the surface density of states (SDOS) at the Fermi level 
when a vortex is situated within a coherence length of the surface of a chiral p-wave 
superconductor. These authors predict a large anomaly in the SDOS depending on the 
sign of the chirality based on a non-self-consistent model for the vortex core which omits 
the internal structure of the vortex cores of chiral superconductors described above. 

Similarly, the authors of Ref. [72] report a dramatic difference in the broadening of 
the energy levels in the vortex core for the m = +1 and m = — 1 vortices Jjj| Their analysis 
is based on an assumed form for the vortex amplitude that does not include the induced 
component of the order parameter in the core with time-reversed chirality. Our results, 
which are based on fully self-consistent solutions for the order parameter, single-particle 
self-energy and spectral density, show only a minor reduction in the spectral width of 
the m = — 1 core state compared to that for the m = +1 core state. This is consistent 
with earlier calculations by Hess et al. J73J [73], who found a small asymmetry in the 
dispersion of the Andreev bound states for self-consistent solutions for singly quantized 
vortices in the chiral ground state. 

The discrete symmetry of the point group is not expected to significantly change the 
structure of the m = — l,p = 1 vortex since the next relevant higher order components 
are m = —l,p = —3, +5 . . .. However, the m = +l,p = +3 vortex allows a higher-order 
subdominant component with p = — 1, which may have a significant amplitude since 
it will be localized relatively close to the vortex center. The corrections to the vortex 
core resulting from higher order harmonics to the basis functions allowed by the discrete 
point group symmetry will reduce the difference in free energy between the two types of 
singly quantized vortices. Takigawa et al. [75] have considered the extreme anisotropic 
limit for the structure of singly quantized vortices and the local DOS by solving the 

jj The m = +1 vortex in the chiral ground state with L° rb = +hz corresponds to the L z = +2 vortex 
of Ref. [72], while the m = — 1 vortex in the same chiral ground state corresponds to their L z = 
vortex. 
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Figure 2. Singly quantized vortices. Shown are the dominant order parameter 
A + (first column) and the subdominant order parameter A_ (second column). The 
color coding is the same as in that of Table 13.11 and the size of the symbols reflects 
the magnitude of the order parameter. In the third column we show both order 
parameter components, A + (black) and A_ (red), as a function of distance along a 
trajectory through the vortex center. In the fourth column the local density of states, 
N(x = 0,y;e), is shown as a function of energy and "impact parameter", y, from the 
the vortex center for y = — 12.5£o to y — 12.5£o and a spacing of Ay — 0.78 ^o- The 
red spectrum showing the zero-energy core state is evaluated at x = 0, y = 0. 

Bogoliubov equations numerically on a lattice with nearest neighbor (NN) hopping as 
well as local NN pairing interactions. 

3.3. Doubly quantized vortices 

As in the case of singly quantized vortices, there are two in-equivalent doubly quantized 
vortices: (i) m — 2, p = 4 and (ii) m = —2, p = as shown in Fig. [31 The latter 
case describes a doubly quantized vortex with a nearly homogeneous core amplitude, 
i.e. a "coreless vortex". Since there is no local phase winding associated with the 
subdominant amplitude, it develops to the asymptotic value for the dominant order 
parameter component as indicated in the third column of Fig. [3j This is similar to a 
domain wall separating the two degenerate order parameter components, except that the 
subdominant component is restricted to the core (see Fig. H]). The large subdominant 
component which "fills in" the core of the doubly quantized vortex substantially lowers 
the free energy for the m = —2, p = vortex by recovering nearly all the lost 
condensation energy from the core of the dominant component. This low core energy is 
reflected in the local density of states shown in the lower panel of column 4 in Fig. 
IU there are very few low-energy, sub-gap excitations associated with pair-breaking 
in the vortex core. The large reduction in the core energy of the doubly quantized 
vortex compared to that of the singly quantized vortex allows for stable lattices of 
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Figure 3. Doubly quantized vortices. The columns and notation are the same as 
that of Fig. [2] Note (i) the absence of a zero-energy state at the vortex center, (ii) 
the small sub-dominant amplitude, A_, for the m = +2, p = +4 vortex, (iii) the large 
sub-dominant amplitude that fills the core of the m — — 2, p = vortex (iv) and the 
correspondingly low density of sub-gap excitations in the vortex core. 



doubly quantized vortices at sufficiently high magnetic fields. The possibility of stable 
lattices of doubly quantized vortices in superconductors with broken chiral symmetry 
was investigated within GL theory. It was shown that doubly quantized vortices can be 
stabilized compared to the singly quantized vortices over a wide range of fields below 
H c2 [76J. The axial symmetry of the current distribution (shown in the lower panel of 
Fig. HJ) leads to isotropic interactions between these doubly quantized vortices, and thus 
a hexagonal lattice structure. 

Broken Axial Symmetry in the Vortex Core The calculation of the structure of the 
doubly quantized vortex with m = +2 and p = +4 shown in Fig. [3] is for a relatively 
small area grid. The result shows an axially symmetric structure in which the local phase 
winding of p = +4 is concentrated as a multiply-quantized axially symmetric vortex. 
However, this structure is unstable. The stable solution for this vortex spontaneously 
breaks axial symmetry, with the p = +4 vortex in the time-reversed phase dissociating 
into four singly quantized vortices arranged around the edge of the time-reversed core 
amplitude as shown in the top-right panel of Fig. [51 The breaking of to C4 symmetry 
is energetically preferred because dissociation of the p = +4 vortex allows the time- 
reversed amplitude in the core to recover to the bulk value and restore most of the lost 
condensation energy of the dominant phase. This reduces the energy splitting between 
the two types of doubly quantized vortices and implies that the region of stability of 
doubly quantized vortices will be comparable for either field orientation, H||L 2 = +hz 
or — H||L Z = +Hz. However, the axial anisotropy of the current distribution of the 
doubly quantized vortex (lower panel of Fig. [5]) is expected to lead to a square vortex 
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Figure 4. Axially symmetric double quantum vortex with m = —2 and p = 0. 
Top-left: the dominant amplitude | A + (R)|/27r /cgT c is strongly suppressed in the core. 
The projection in the plane is a phase plot of sgn(?le [A+(R) exp(im</>)]). Top-right: 
the amplitude, | A_ (R)|/27rfcsT c , develops to its equilibrium value and fills the vortex 
core, while the phase plot of sgn(^He [A_(R) exp(ip<p)]) shows a radial tt phase change 
indicative of a domain wall separating the time-reversed phase in the core. Bottom: the 
domain wall separates an axially symmetric flow in the core that is counter-circulating 
relative to the current outside the core. 

lattice of doubly quantized vortices for H||L 2 = +hz, compared to a triangular lattice 
for — H||L a = +hz. 

Finally, let us consider the effects on vortex structure from additional order 
parameter components which are allowed by discrete point symmetry of Fermi surface. 
For the axially symmetric double quantum vortex with m = —2 and p = the next 
order harmonics are those with p = ±4. The magnitudes will be very small because both 
amplitudes, A + and A_, are large and there is no energetic advantage to inducing these 
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Figure 5. Spontaneously broken axial symmetry in the core of the double quantum 
vortex with m = +2 and p = +4. Top-right: The p = +4 vortex in the A_ dissociates 
into four singly quantized vortices situated as satellites at the corners of a square on the 
core boundary (note also the phase plot in the projection), allowing the amplitude | A_ | 
to grow substantially in the core. The satellite vortices of A_ lead to 4-fold anisotropy 
in the dominant amplitude |A+| (top-left panel), as well as a current distribution 
(bottom panel) which exhibits the four-fold anisotropy induced by the four p = +1 
core vortices. Note that there are no counter-circulating currents once axial symmetry 
is broken. 

components. Thus, we can safely neglect these higher harmonics. The same reasoning 
applies for the vortex with m = +2 and p = +4. Both the dominant and core amplitudes 
are large and suppress the higher harmonics, even the higher harmonic with n = —4 
(i.e. p = 0). This would not be the case if the axial symmetry were not spontaneously 
broken. 
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3-4- Structure of inhomogeneities in the Meissner phase 

In addition to novel vortex states, broken chirality leads to defect structures that carry 
current, but have no global phase winding. The example of a cylindrical normal metallic 
inclusion is included in Table 13.11 This defect has no global phase winding, i.e. m = 0, 
but necessarily has a local phase winding of p = +2 for the induced order parameter. 
Thus, in contrast to vortices that are stabilized by a field, this defect can exist in the 
Meissner phase. Such defects can be engineered, or as in the case of Sr 2 Ru04, Ru ions 
may precipitate and form metallic inclusions embedded in the superconducting host 
material. 

The ground state of a chiral, p-wave superconductor will be either the p x + ip y 
state, or the — ip y state, although the actual zero-field superconducting phase may 
include a number of domains of different time-reversed states separated by domain 
walls. Here we consider a normal metallic inclusion of radius r incl = 0.78£o embedded 
in the homogeneous p x + ip y ground state, or at least several coherence lengths away 
from a domain wall. We assume a simple model for a metallic inclusion in which 
electrons and holes are perfectly transmitted across the (NS) interface between the 
inclusion and the superconducting material when the latter is in its normal state. We 
also assume there is no pairing interaction within the metallic inclusion. Thus, a mean 
field order parameter does not develop inside the inclusion, however pairing correlations 
and associated spontaneous currents do develop within the inclusion as a result of the 
proximity effect and the high transmission of the NS interface. More detailed models, 
including specific models of Ru inclusions, which incorporate the finite reflectivity of 
the NS interface, as well as pairing attraction within the metallic inclusion can be 
implemented, but are outside the scope of this article. 



A, A_ 




Figure 6. The structure of a normal metallic inclusion embedded in the (p x + ip y ) 
phase. The inclusion is located in the center in column 1 and has a radius r incl = 0.78£o- 
The columns and notation are the same as that of Fig. [2] Note (i) the phase winding 
of the time-reversed order parameter (column 2) and (ii) the shallow bound states 
induced by the inhomogeneous order parameter near the boundary of the metallic 
inclusion (column 4). 

Shown in column 2 of Fig. [6] is the time-reversed component of the order parameter, 
with a phase winding of Air, that nucleates near the metallic inclusion. Broken time- 
reversal symmetry is revealed by the appearance of spontaneous supercurrents that 
flow in and around the metallic inclusion as shown in the top-right panel of Fig. [7J 
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Figure 7. Spontaneous local magnetic field (top-left) and current density (top- 
right) around a mesoscopic metallic inclusion with phase windings of m = 0,p = —2. 
The change in sign of the induced field is indicated by change in color from blue to 
red. For comparison we show the current densities for a singly quantized vortex with 
m = — l,p = 1 (bottom- left), and the doubly quantized vortex with m = — 2,p = 
(bottom- right). 

These currents generate a local magnetic field, b = b(x,y)z, which oscillates in sign 
as one moves radially away from the defect on the scale of the coherence length, £o, 
shown in top-left panel of Fig. [7J The spatial average of the field is zero, but the 
locally varying field is non-zero and should be observable, e.g. as broadening of the 
/iSR linewidth below T c . We also note that the current and field induced by a non- 
magnetic point impurity in a chiral p-wave superconductor [77J is similar to that of a 
mesoscopic metallic inclusion, however the origin and interpretation of the structure 
appears different. Complex current and field distributions may also be generated in 
unconventional, but non-chiral, superconductors by magnetic and spin-orbit scattering 
impurities [19] . 

Finally we note that the leading higher order harmonic corrections resulting from 
discrete C4 lattice symmetry are expected to be weak for a small cylindrically symmetric 
defect. The most important higher harmonic will be a subdominant component with 
p = —2, but it is weaker than the leading p-wave harmonic with p = +2. All other 
harmonics can be neglected for small inclusions because the higher winding numbers will 
force the induced amplitude to be further from the vortex core and thus more effectively 
suppressed by the dominant m = order parameter. A possible exception may occur 
for large inclusions, r incl ^> £ - m this limit the harmonic with p = 6 may be more 
important than the p = —2 component. 

4. Conclusions 

In conclusion we have analyzed and carried out self-consistent calculations of the 
equilibrium structure and excitation spectrum of vortices in layered, spin-triplet, p-wave 
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superconductors and superfluid 3 He films with spontaneously broken chiral symmetry 
in the ground state. We show that the cores of vortices contain the time-reversed 
phase with a local phase winding that differs by ±2 x 2n depending on the chirality, 
i.e. L° rb = ±h, of the ground state. This implies that vortices with equal but opposite 
phase winding are inequivalent, and thus the lower critical field for nucleation of vortices 
in layered, chiral p-wave superconductors depends on the relative orientation of the 
field and the orbital angular momentum of the ground-state pairs. Axially symmetric 
doubly quantized vortices with zero phase winding in the core are predicted. A lattice 
of these vortices is energetically favored at sufficiently high fields for — H||L 2 = Hz. For 
the opposite phase winding, i.e. H||L 2 = +z, the doubly quantized vortex contains 
four circulation quanta for the core amplitude. This vortex spontaneously breaks axial 
symmetry by dissociation of the core vorticity into four singly quantized vortices in 
the time-reversed order parameter induced in the core. A lattice of these vortices 
is also expected to be energetically favored at sufficiently high fields with a square 
lattice structure. In addition to vortex states, inhomogeneities in the distribution of 
non-magnetic impurities, or metallic inclusions embedded in the ground state lead to 
spontaneous supercurrents and magnetic fields induced in the Meissner phase, which are 
localized near the inhomogeneity. These currents are carried by bound electronic states 
associated with an induced order parameter with phase winding of 4tt. The magnetic 
field and current distribution vary on the scale of the coherence length and are expected 
to be observable by sufficiently small local magnetic probes, or as a broadening of the 
/iSR linewidth below T c . 
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